A Model for Dynamic Simulation and Analysis of Tether Momentum Exchange 


Stephen Canfield*, David Johnson 
Department of Mechanical Engineering 
Tennessee Technological University 
Cookeville, TN 38505 


Kirk Sorensen, Ken Welzyn 
TD40,TD15 

NASA Marshall Space Flight Center 
Huntsville AL 35812 


Abstract: 

Momentum-exchange/electrodynamic reboost (MXER) tether systems may enable 
high-energy missions to the Moon, Mars and beyond by serving as an "upper stage in 
space". Existing rockets that use an MXER tether station could double their capability to 
launch communications satellites and help improve US competitiveness. A MXER tether 
station would boost spacecraft from low Earth orbit to a high-energy orbit quickly, like a 
high-thrust rocket. Then, using the same principles that make an electric motor work, it 
would slowly rebuild its orbital momentum by pushing against the Earth's magnetic field- 
-without using any propellant. One of the significant challenges in developing a 
momentum-exchange/electrodynamic reboost tether systems is in the analysis and design 
of the capture mechanism and its effects on the overall dynamics of the system. This 
paper will present a model for a momentum-exchange tether system that can simulate and 
evaluate the performance and requirements of such a system. 


1. Introduction: 

Many uses for space tethers have been proposed and discussed in the literature (see for example 
Beletsky and Levin [1], Carroll [2]). One area of proposed use is in creating an in-space payload 
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transportation system. For example, Carroll [2] presented an early tether momentum transfer system that 
functioned by changing the tether length in resonance with the in-plane libration frequency to accelerate a 
payload. Hoyt, [3] and Forward [4] demonstrated a concept for performing momentum exchange 
between a spinning tether and a payload to transfer the payload into a lunar orbit, with capture and release 
of the payload performed at perigee of an elliptical tether orbit. Bangham, Lorenzini and Vestal [5] 
analyzed both single and multiple spinning-tether systems for boosting payloads from low-earth to 
geostationary orbits, with the multi-tether system offering the lowest mass. Hoyt et al. [6] and Hoyt [7] 
demonstrated a single tether system for boosting payloads from a low Earth orbit to a geostationary orbit. 
This system was designed to also provide transfer of payloads to limar orbits and future transfer to the 
Moon and Mars. Sorenesen [8] considered the results of the Tether transporation study (Bangham, 
Lorenzini and Vestal [5]) and concluded that a non-equatorial, multi-tether system is nonviable due to 
orbital regression. Sorensen also discussed the use of electrodynamic reboost for the tether facility, 
described the tether/payload rendezvous, and presented a passive capture mechanism for the rendezvous 
process. 

One important issue in the early consideration of a tether momentum-exchange facility is the dynamic 
behavior of the system. Several tether dynamic models have been presented in the literature, with a 
number of commercial or applied tether model analyses developed. Misra and Modi [9], offer a thorough 
review of the modeling approaches for electrodynamic tethers while continued work includes lumped- 
bead formulations (Keshmiri and Misra [10]) hinged-rod models (Puig-Suari et al [11], Biswell et al 
[12]), analytic models (Beatty and Haddow, [13]) and numerical techniques for simulating system partial 
differential equations (Yokota, Bekele, and Steigmann [14]). These models primarily consider single or 
multi-chain electrodynamic tethers during deployment, flight and retrieval from a larger vehiele. In the 
application of tether momentum exchange, a simulation of the nonlinear dynamic behavior of a free, 
spinning tether system in orbit and during the momentum exchange process is of interest and the subject 
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of this paper. This paper will investigate the dynamics of the momentum exchange process based on a 
spinning, orbiting tether facility and a passive capture mechanism as suggested by Sorensen [8]. The 
paper will develop a model for this tether system based on a finite element approach and perform 
subsequent simulation and analysis of the dynamics of the momentum exchange (MX) procedure. This 
approach is based in part on the techniques of finite element analysis of dynamic systems demonstrated 
for example by Meirovitch [15] and Nath and Ghosh [16]. The model will be general to allow a variety 
of tether and capture element types to be included in the future, allowing the simulation to be used as a 
design-guide in determining the parameters in a capture device. The model will be employed to simulate 
the rendezvous process with the effects on the tether system examined and considered. 

2, A Model for Tether Simulation 

The dynamics of a flexible tether in space flight are evaluated using a finite element model 
formulation of the continuous flexible tether system that will include expressions for the coriolis, 
tangential and normal components of elastic accelerations. The objective of this dynamic model is to 
observe the time response of the flexible system during tether flight and through the momentum exchange 
process. This finite element dynamic model will represent the tether system as a series of cable elements, 
each undergoing a combination of rigid body motion and axial elastic deformation and allowing free 
rotation between adjacent cable elements. The finite element method will be applied to this system as a 
means of representing the total tether motion as the motion of a continuous assembly of finite elements, 
with each element representing a part of the tether system, and requiring continuity and compatibility 
between connecting elements. 

2.1 Finite Element Model 

The finite element approach provides an approximate means to express the displacement of any point 
within a continuous element in terms of the displacement of a finite number of nodes multiplied by a 
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suitable set of interpolation functions (Meirovitch [15]). The interpolation functions are typically selected 
as a polynomial of order n, generally of low-degree, £ind are the same for all elements of a given type. 
Use of these interpolation functions to describe element motion can then result in approximate solutions 
for the governing partial differential equations for each element. 

Consider the case of a single tether element in an equatorial orbit around the earth, shown in Fig. 1. 
This element represents a small portion of the entire tether system divided into a series of similar 
elements. The equations of motion will be derived for this representative tether element using the finite 
element method approximations and Lagrange’s equations, and will provide the basis to assemble a 
system of similar elements and thus a complete tether system. 



Figure 1: Tether element in equatorial Earth orbit 

The location of particle dx on the tether element shown in fig. 1 is given by Xp as: 

Xp (x,/) = X, (/)-h (x + u(x,t)) Eq. 1 
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with u(x,t) representing the elastic axial displacement of the element at point x along the element in time, 
both described in the element frame, and Rekm a suitable rotation operator expressing the displacements in 
a uniform frame. Using the selected interpolation functions, the axial displacement can be written as 

w(x,/) = l(xf u(/) Eq. 2 

with l(x)^ = [i, (x) Lj{x) ...]^ the interpolation functions and u(/)=[m|(/) ^ 2 (f) ‘"F nodal 
displacements in time. The particle position can be defined as a combination of these interpolation 
functions and nodal positions as, 

Xp 0 = X, (0 + R./.M (x + l(x7u(/)). Eq. 3 

A set of generalized coordinates will be defined to describe the motion of the tether element. This set will 
consist of the reference position, Xr, the orientation of the element defined in a rotation operator, R, and 
the nodal displacements, u. This vector of generalized coordinates, q, becomes 

q = R|i Eq. 4 

with R containing element rotation and p. defined as 

= 1 M, Hjf. Eq.5 

The position of any point along the element is then given as a combination of the interpolation functions, 
cast into matrix form, and the generalized coordinates, 

x^=L^Rp. Eq. 6 

Note that contains the spatial dependencies along the tether length while R and p are functions of 
time. To proceed with Lagrange’s equations, the Lagrangian must be formed as the combination of 
kinetic energy and potential energy in the tether element. The kinetic energy is first defined over the 
element. 
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with / the tether element length AfxJ the element cross section area and pfxj the element density. The 
velocity of the tether element at location x is 

i^=L^R^ + L"RH Eq. 8 

with R providing the angular velocity of the tether element frame. This will be recast to combine the 
generalized coordinate velocities into a single vector, |i' as 

Xp=L'"R>' Eq.9 

where R’ is a general operator that gives velocity orientation. 

This paper will consider tether motion in an equatorial plane. Figure 2 shows the generalized 
coordinates for this planar case, in which |i' , and R’ for element 1 become 

1 0 0 0 0 0 0 

0 1 0 0 0 0 0 

0 0 0 0 0 0 

^riOxOZ, OZ^Ol^OOci? 0 0 0 0 

[o 1 0 X 0 I, 0 ij’ 0 0 0 -S0 0 C0 0 

0 0 0 c<9 0 0 

0 0 0 0 -s0 0 c0 

0 0 0 0 c0 0 s0 
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Figure 2: Generalized coordinates of a tether element 


Using the expression of given in Eq. 9 , the kinetic energy can be rewritten as, 

7 = ^1 p(x).4(x)(l^R' |i') (ix = ^ R'^ | p(x)^(x)LL^atcjR'p' Eq. 1 2 

or 

7 = ^|i'^R’^iii^R'|i' Eq. 13 

with me the consistent element mass matrix determined from the interpolation matrix and element area 
and density as 

= |/?(x)y4(x)LL^cic . Eq. 14 

This element mass matrix me is shown for reference in the appendix. In a manner similar to evaluating 
the kinetic energy, the potential energy of the tether element is evaluated. In this derivation, the potential 
energy will consider strain energy in the element only, while gravitational potential along with other 
external loads will be included in the generalized force. The potential energy in the element is given as 
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Eq. 15 


Using the expression for Xp given in Eq. 6, 

V = - dx = -n'"R''k,Rn Eq. 16 

I 5x j 2 

with A and E the element area and modulus, and k* the element stiffness matrix evaluated as, 

k, = lAE(x)^^dx. Eq,17 

The element stiffness matrix k« is given for reference in the appendix. The generalized force is 
constructed next to represent external forces located at the nodes. This is derived from an expression of 
the virtual work resulting from the gravitational loads and all nonconservative forces distributed over the 
element, 

SJV - = |f(x,t)-L^j(Rn)c6c = F^(Rfi) Eq. 18 

with f(x,t) the external force distribution over the element and F the external nodal force given as: 

F^ = |f(jc,/)L’’o!x. Eq. 19 

Partial derivatives of the Lagrangian, 3 = T-V , with respect to the generalized coordinates and their 
time derivatives are evaluated next. In general form, the elements of Lagrange’s equations are considered 


as. 



R’^ m,R>' 


Eq. 20 



dt 








,. R '^m.R>’+ 






Eq.21 
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Eq. 22 

and the generalized force associated with qt is derived from the virtual work as. 


Eq. 23 


From these partial and time derivatives, the equations of motion for the element are assembled for each 
generalized coordinate, /, 


1-^- =a 


- R"- m.R'pi’-|i'' m,R>'+«^^^ R’'k,R|. + R"®^^ k.R|i = Q, 

for 1=1 ... m, to result in m equations of motion with m the number of generalized coordinates. 


Eq. 24 


Eq. 25 


After expanding the system of equations in Eq. 25 for each generalized coordinate and combining 
accelerations, the equations of motion take on the general form 

Meii'+DXg.^V+KX^V = Q* Eq. 26 


with 


M, 




R" m R’ 



Eq. 27 
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Eq. 29 



Eq. 30 


An expansion of each of these terms is listed in the appendix. Thus, Eq. 26 provides the equations of 
motion for the tether element of Fig. 2. It now remains to perform an assembly process by which the 
equations for all the elements forming the tether system are compiled. 


2.2 Element Assembly: 

The assembly procedure allows the entire tether system to be represented as a collection of elements, 
with geometric compatibility maintained at all element nodes. Let fig. 3 represent a portion of the tether 
system modeled as a collection of tether elements. 
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Figure 3: Collection of elements forming the tether system 


Each element is attached to the end node of the previous element and maintains continuity through 
nodal displacement. As additional elements are added to the tether system, additional generalized 
coordinates appear. For example, in this case of a planar tether system, two new generalized coordinates 
occur with each additional element. Thus, the vector of generalized coordinates for a tether with n 
elements is, 

^ = »1 «2 ^2 T ' Eq. 3 1 

For a system with n elements, there are m = 3+2« generalized coordinates. 

The equations of motion for the entire tether system are constructed through assembly of the individual 
element equations of motion, while taking into account continuity at the nodal positions. The assembly 
procedme can be observed by forming the total kinetic and potential energy and virtual work of the tether 
system as a sum of the energies and virtual work contributed by each element. 


^ 1 = 1 


1 ^ 


Ei 4 £,(^ 


Ke,Xer=-Z2, 
^ /=! 


'ai/Rj*,^’ 

Sc 


dx 


^>f',.„=tlt,(x,iyL'S(R,ll,)jx 


which leads to Lagrange’s equations for the system, 


— 7 ’ _y ^^telher 

"^tether ^ tether ^ tether ’ 


Eq. 32 
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To make the individual elements compatible in the summation processes, a uniform set of descriptive 
coordinates must be developed for each element, similar to Eq. 6 for the first element. Consider an 
internal element, call it the t'* element, located along the tether (refer for example to Fig. 3). The position 
of a point along the i'* tether can then be given using the same form as equation 6, 

Xp, = L^R,|i, Eq. 34 


with the displacement vector, pj given as: 


F/ = 


+^cos 1 u, M,.,, 


Eq. 35 


L M M J 

and R| the rotation operator for element i. Thus, the position of a point on the tether element is a 
function of the tether reference vector as well as all preceding nodal displacements and element 
orientations. In a similar manner, the velocity of a point on the /** tether element is generalized to 
become. 


Xp,=LX,iV 


Eq. 36 


with R’l a general operator that gives velocity orientation for element / and 
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i^+SkcoX'^jX'y+Vil+^C^yK.] “< “- 

/-i y-i 


Eq. 37 


Using these general expressions for the position and velocity of reference positions for each element, the 
equations of motion can be derived for the entire system through the summation process demonstrated 
above, resulting in a general form for the equations of motion for the entire tether system, 

Mq + D{q, q)i{ + K{q)q = Q Eq. 38 

where q , q and q represent the generalized coordinates associated with the entire tether system, and given 


in the general vector form as, q = 6^ u^ &2 "a •r and the matrices, M, D, and 

K and vector Q are coefficients of the generalized coordinates and their derivatives. These are found in 
general through expanding the elements of Lagrange’s equations for each element, shown for example for 
element 1 in Eqs. 20-23, and then combining terms to result in the necessary coefficient matrices of the 
generalized coordinates. The procedure is demonstrated in the appendix, with the M, D, and K matrices 
and the vector Q expanded for a two element system. Note that since element rotation is decoupled from 
the stiffness in the case of geometric linearity, the ^associated terms in the generalized vector multiplied 
by the stiffness matrix will be set to zero. With the general equation of motion generated for the tether 
system, a solution approach used in tether simulation is presented next. 


2.3 Solution of Equations of Motions 

Eq. 38 provides a system of second order differential equations of motion of the tether system defining 
the response of the generalized coordinates in time. Combining the shape fimctions and the nodal 
displacements in time, the approximated response of the continuous tether system is observed. This 
behavior of the rotating tether system in orbit is the point of interest in this work, and is solved using a 
direct integration routine. A Runge-Kutta integration procedure with adaptive time step is selected and 
employed for dynamic tether analysis. 
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2.4 Element types used in the Tether MX model 


2.4.1 Tether Element: 

An element with a single longitudinal non-rigid-body degree-of-freedom is first employed in the planar 
tether model. Elemental longitudinal stiffness and damping are associated with this non-rigid-body 
degree-of-freedom. Each tether element is modeled with two nodes representing a linear shape function. 
Rigid body motions of the tether element include rotation and translations in the orbital plane. Due to the 
inability of the tether element to support compressive loads, a significant nonlinearity occurs when the 
tether goes slack. Several models have been proposed to accommodate this nonlinearity, the simplest 
approach is employed in this model by defining the element stiffness to go to zero when the element goes 
slack. 


^i—hhn 

[O. /, < 


Eq. 39 


Note that in the spinning tether configuration, the tether elements maintain positive tension at all times 
due to centripetal acceleration forces and slack conditions are generally identified as unstable conditions 
in the tether model. 


2.4.2 Ballast and Capture Element: 

The ballast or facility element represents a rigid body, ballast mass located at one end of the tether. 
This element will be composed of an infinitesimally small member with mass localized at a common 
node. Note that since the length approached zero, the stiffness approaches infinity, this is represented in 
the global matrix with a single node. The mass of this element is also associated with a single node, and 
the force effects of this element are derived for the single nodal position. 

2.4 Tether Momentum Exchange 
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Momentum exchange occurs through the process of capture, cany and release of the payload. Initial 
capture of the payload is modeled as an inelastic impact process between the last node of the capture 
element and a rigid lumped-mass payload. The rendezvous process is included in the tether model in the 
following manner. The tether model is simulated through integration forward in time starting from a set 
of initial conditions to result in a future rendezvous with the payload. At the point of rendezvous, time tr, 
algebraic evaluation of the new velocities due to impact between the point mass payload and the mass 
localized at the end of the tether is performed. This results in a new set of boundary and state conditions 
that are used to update the tether model and continue the time integration of the tether model. Payload 
release involves again an instantaneous change in the boundary conditions of the tether model. Note that 
in general, the rendezvous process is considered to occur with the end of the capture device and the 
payload having the same position and velocity, resulting in no net effect from impact. Error in either 
position or velocity of the tether or payload would add additional d)mamic effects to the system. 

The remainder of this paper reports on some of the results from the simulation of the momentum 
exchange process. 

3.0 Analysis of Tether Capture and Release of Payload: 

Using a simple FEA tether model, the capture process is simulated as follows. The tether simulation is 
started and proceeds for a period of time sufficient to allow initial condition transients to die out. 
Similarly, the payload is released and propagated in time in its proper orbit. Initial conditions of both the 
tether and payload are selected to result in a rendezvous precise in space and time. The momentum 
process begins when the payload enters the capture region of the tether and continues through the cany 
and release process. For the general case, the capture oceurs when the tether is at perigee and the tether is 
aligned with the local vertical, the capture end pointing toward earth. The tether than carries the payload 
through a predetermined angle and releases the payload, with both the payload and tether assuming a new 
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trajectory. The simulation process records the data from the tether model including tether tension, 
rotation rate, center of mass, and velocity. These results provide essential information in predicting 
momentum exchange performance under ideal conditions, setting the bounds on required tether 
propagation accuracy, and in providing design guidelines for design of the capture device. Results from 
analysis in each of these areas follow. These results are based on tether and payload sizes published by 
Hoyt [7], as shown in Table I. 


Table I: Tether and Payload Parameters 


Constants 


Tether Facility Orbit 



3.986E+14 



Payload Mass (kg) 

1000. 2500 

Apogee (m) 



250 

Time for one orbit (s) 

11,176.25 

Ballast Mass (kg) 

11000 

SMA 

10,804,136.00 


80000 

e 

0.371987172 


1.9635E-05 

Perigee Velocity (m/s) 

8977.690593 


970 


4109.444293 



I'l 111 llil 111 1 

-0.013881192 

Orbit Radius (m) 

6694275,739 



Vp (m/s) 

7716.44186 




5450.876936 








3.1 Results From Tether Momentum Exchange Simulation 


The momentum exchange analysis will demonstrate simulation of tether and payload rendezvous, 
capture, carry and release, under ideal conditions. Here, ideal conditions imply that there is no error in 
position or velocity of either the tether or payload. At the point of rendezvous, both the capture device 
and payload have the same position and velocity. The tether produces and instantaneous application of 
acceleration on the payload in the direction of the tether, resulting in increase in tether tension with much 
smaller effects on transverse vibrations. The results from this ideal capture are shown in figs. 4-8. 
Figures 4-6 show the tether orbit radius before and after rendezvous, with Figs. 5,6 zoomed to show 
effects of capture and release. Figure 7 shows the payload orbit before and after rendezvous while fig. 8 
gives a time history of the tether tension. Finally, Table II below summarizes the results of the orbits of 
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the tether facility and payload before and after the momentum exchange process. Table II shows the 
results for both a 1000 kg and a 2500 kg payload. 



Fig. 4: Orbit Radius of Tether Fig. 5: Orbit Radius of Tether 



Fig. 6: Orbit Radius of Tether Fig. 7: Orbit Radius of Payload 
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Tension in Element 2 


X 10* 



'lOOO 2000 3000 4000 5000 6000 7000 6000 9000 10000 


Time 


Fig. 8: Tension in Tether elements 
Table II: Summary of Tether and Payload Orbits 



Payload 

Tether Facility (1000 kg Payload) 

Tether Facility (2500 kg Payload) 

Pre Catch 

Post Release 

Pre Catch 

Post Release 


Post Release 

Perigee 

6693863 

6789026 

6785077 


6785077 

6763012 


315727 

410890 

406941 


406941 

384876 



97999789 



14823136 

BKSBHI 

Apogee Altitude 

SESSHHHI 

91621653 

8445000 


8445000 


Period 


119354 


9653 

11176 

8357 

Seml-maior Axis 

6693863 

52394407 

10804106 

9798655 

10804106 



0 

0.8704 

0.3720 

0.3095 

0.3720 

0.2402 


7717 

10479 


8783 

8978 

8549 

Apogee Velocity 

7717 

726 


4632 

4109 

5238 

U) (deg) 



^^^^||||||||||[||||^ 

2.75 

0 

3.63 


From the error free momentum exchange case several results are noted. First, general effects on the 
tether orbit correspond to those published in the current literature (Hoyt [7], Bangham et al [5]). 
However, the tether orbit is also affected by elastic material change during the capture and release 
process. The tension undergoes a significant increase during both capture and release. The tether tension 
is fairly stable before and after capture with small oscillations in amplitude due to tether rotation in the 
gravitational field. At rendezvous, the tether experiences a significant spike in tether tension. Release 
again causes a step input to the tether system, resulting in a damped transient response returning to the 
free-tether steady-state conditions. The response of the tether during the momentum exchange process is 
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largely dependent on the stiffness and damping characteristics of both the tether and capture elements. 
An additional effect is the resulting tether orbit. As predicted, the tether facility loses altitude, primarily 
at apogee. The orbit change in the tether facility for both the 1000 kg and 2500 kg payload is significant, 
but gives the tether sufficient opportunity to reboost for future momentum exchange. A final 
characteristic of the system is noted. The line of apisides of the tether is shifted at MX due a finite time 
between capture and release and the elastic behavior of the tether system. 

4.0 Summary 

The momentum exchange process has been modeled using an FEA approach, with the basic processes 
verified against current literature. The FEA model shown is based on the simplest single degree of 
freedom tether element and accoxmts for dynamic effects due to the rigid body motion and elastic 
deformation of the element. This model can be extended to multiple tether elements and allows for 
specification of various element types. The paper then presents an analysis of the dynamic effects of the 
momentum exchange process on a tether system and payload. This analysis demonstrates relative 
stability in the process and confirms many of the previously reported results. It also demonstrates a 
number of dynamic effects in the process such as longitudinal tether stretching and a shift in the line of 
apsides. Finally, the paper demonstrates the ability of the tether model to boost a variety of payloads 
given a sufficient time for reboost. 
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Appendix: 


The general displacement of a particle along element i is given as a combination of the interpolation 
functions in the form. 




with L, R|, and pi given as, 


= 


+ ^cos JC^+£sin((0-X/;+«,„) 1 «, 




y=i 


Eq. A1 




1 0 :c 0 ^ 0 ^ 

0 1 0 JC 0 0 


1 0 0 0 0 

0 10 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 c^, 0 

0 0 0 5^, 0 

0 0 0 0 c6», 

0 0 0 0 s0, 


Eq. A2 


with s&i, cOi = sin(6j), cos(^) respectively. Similarly, the particle velocity will be cast in the general form. 




Eq. A3 


with 


l»r 


U„ + 


y=i y=i 
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Based on Eqs. 14 and 17, the local mass and stiffness matrices for each element are given as. 
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and the local force vector is given as. 


Fe = 




fy {fJ 


fyl 


fj \fj \fyl 
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with fx,fy the force components acting on element i. 

The elemental mass, stiffiiess matrices and velocity product matrices, M^, Kei, Dei resulting from the 
components of Lagranges equations given in Eqs. 20-23 are found as, 


M„. = 


/ 




I T 


dq, 
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^diiry ]_dii\y 


VL 


dt 
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|k„R, 


jy 


Eq. A6 


while the elemental generalized force vector is, 


Q„ = 



Eq. A7 


The system mass, stiffness and velocity product matrices are formed as combinations of the elemental 
matrices in the fashion demonstrated in the following equations. These matrices contain the coefficients 
of the generalized coordinate vector and its derivatives in the form 

g = [x„ 0 , M, M2 02 «3 -J. 


The global meiss matrix, M is generated as. 
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/=1 

M(:,2)=XM„.(:,2) 

M(:,3) = M,, (:,3)+ «,M„ (;,4)+ (:,5)+ JCA + M 2 X- c,M„.(:,2)) 

1=2 

M(:,4) = M,.(:,6) 

M(:,5) = M„ (:,7)+ M.,(:,6)+ X(c,M„(:,l)+:?,M,X:,2)) 

i=2 

M(:,6) = (:,3)+ (:.4)+ «3M,2(:,5)+ ^(/j + MjX" ^ 2 ^^ (:,1)+ (;,2)) 

1=3 

The global velocity product matrix D is generated in a similar form as, 

d(:,1)=ZD,(:.1) 

i=\ 

D(:,2)= Jd„.(:,2) 

i = I 

D(:,3) = D„ (;,3)+ m.D,, (:,4)+ M3D„(:,5)+ ^(/, + « j5,D^(:,l)+c.D„(:,2)) 

1=2 

+ Z(m„(:,1)(-^iC,(/, +Uj)-2«,5,)+M^(:,2X- 021^1 (^i +«2)+2MjC,)) 

1=2 

+ M,M„(;,4)+ii2M^i(:,5) 

D(:,4)=D,,(;,6) 

D(:,5)=D.,(;,7)+D.,(:,6)+XM„(:,1)+^,D„(:,2)) 

/=2 

D(:,6)= Dj 2(->^)'*’ “2®e2(^^)'*' ]^(^2 X‘*2®ri(->0'*' ^2®« (->2)) 

1=3 

+“3)-2mj^2)+M„(:,2)(-0j5j(/2 +M3)+2M3cJ) 

(=3 

+ MjM,j(:,4)+U3M,j(:,5) 


Eq. A8 


Eq. A9 
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The global stiffness matrix, K is generated as, 


K(:,1)=2;K„(:,1) 


/=! 


K(:,2)= tK„(:,2) 


/==! 


K(;,3)=:0 

K(:,4)=K„(:,4) 

K(:,5) = K,X:,5)+K,,(:,4)+X(K„(U)c,+K^(:^^^^ 

K(;,6) = 0 


i=2 


Eq. AlO 


and the generalized vector of forces, 


n n 

q=EQ. = Z 


i=\ 1=1 


F„. R.. 




«i — , /Qq 




0R, 






Eq. All 


The system mass, stiffness and velocity product matrices and the total generalized force vector are 
developed as an example for a simple two-element system. This system has three nodes, two elements, 
and seven generalized coordinates. The results of this expansion are given in the following equations. 


M: 


Eq. A12 


M(l,:) = [2/p 0 ^Vp(^ + “i + 3 M 2 ) 2<^i^P T^pC^i+^Cj) ^Cj/p] 

M(2,i)=[ 0 2/p yCj/p(3/ -1- Mj + 3 ^ 2 ) ' 2 ‘^l^P '2 ^pC^I 4" 3 ^ 2 ) Y ^ 2 ^P(^ 4" 1^2 4" W 3 ) '2‘^2^p] 


M(3,:) = 
M(4,:) = 


f^sjp{l + u^+3>u^) 4.c,/p(3/ + M, + 3 U 2 ) w*33 m 

jc,/p |5,/p 4 -c, 5,/^P i/p i/p 0 0] 


m 


37 
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M( 5 ,:)=[i/p( 3 c,+c,) i Ipi^s^+s^) W53 \lp |//?(7 + 9 c, Cj +95,53) Wjj |/p (7 + 3 c, C2 +35,53)] 

M(6,i)= — -j53/p(~/ + M2 +Wj) yC2/p(/ + M3 +Mj) ^53 0 -JC2S2I P ^66 T*'2‘^2^ p] 

M(7,:)=[^C2/p j53/p |(/ + M2)/p(-5,C2 +c, 52) 0 |/p \C2S2l^9 j/p] 


with, 

WI33 =j/p(4 /^ -/m, + 2 lu^cf +4/M2 + 4 lU 2 cf +Mj* +M,«2 +4M2) 

MI35 =i/p(4c,/5, -95, /C2 -95, M2C2 +9c, /52 +9c, M252) 

/W 36 =i(/+»,Vp(- 5 , 53 / + 5,52M2 + 5,52M3 + C,C 3 / + C,C 2 M 2 + C,C 2 Mj) 

^37 =2(^ + "2Vp(-^1<^2 +^1*^2) 

/Mj 3 =j/p(4c,5,/-35,C3/ + 3c,52/-35,C2M2 +3c,52M2) 

"*56 = 6^P(“^^1‘^2^ -252C3/-35,C2/ + 3 c, 52M2 -35,C2M2 +3c, 531/3 -35,C2t/3) 

"*63 =y(/ + M2Vp(“'^1‘*'2^ + ‘^1'^2“2 + ■S'i-S' 2“3 + C,C3/ + C,C2M3 +C,C2M3) 

"*66 —3^9^^ ~ ^“2 + 2cllU2 - 2/M3 + 4C2/M3 + Mj + M2M3 + Mj ) 

D: Eq. A13 


D(l,:) = 

D(2,:) = 
D(3,:) = 
D(4,:) = 
D(5,:) = 
D(6,:)= 

D(7,:) = 


[0 0 y/p( 20,C,/ + 20, M3C, +51/35, +M,5|) 0 0 -j/p52(2M3 +1/3) 0] 

0 0 ^/p(20,5,/ + 20,M35, -51/2C, -M,C|) 0 0 7/pC2(l/2 +M3) 0] 

j0/|C|/p(/ — M, — M2 ) y0/,5,/p(/ + M, + M2) </33 //34 //35 6/35 oj 

^/p0,5, y/p0,C, |/p0, (-/ + 2/c,^ +2 m, +M 2) 0000 ] 

f/p(0,5, +0252) ^/p(0,c, +03C3) J53 0 0 0] 

:^02C2/p(-/ + M2 +M3) f0252/p(/ + M2 +M3) 0 d^^ d^ d„] 

f 0353/p I03C3/P d^^2 0 0 j03/V(-l + 2C3^)+i02/p(«2+2M3) 0] 


with. 


d^y =^/p(4M|0,5, +80,M25,c,/-1Om3/-4c,^/m3 -I 4M2M2 +m,/-2m,/c,^ -2m,m, -MjMj -M3«,) 
J 34 = j0/,/p(2cj*/ - / - 2 m, - M 2 ) 

^35 = -11 -u^- IU2 ) 

6/35 = y/p(4M2/5,52 +M2M25,52 + M2/C,C2 + M2M2C,C2 + 1/3/5,53 +M3M25,52) 


^53 


-4/0,C,^ +80,/-M,0, +40, M2 + 35,5262/ -3 c,C 202/ + 35,5202M2 -3 c,C 202M2 

+ 3c,C20,/ + 3 c,C 20,M2 + 65,C2l/2 + 35,520,/ + 35,520,M2 - 6c,52M2 
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Jj6 =j/p(02^ + 262/^2 +2 i/ 202 +«302 -3li2^iC2 + - 3 u ^ C ^ S ^ +3M3^,C2) 

"^63 = "2 (" 2 ^1^2® 2^p(“ / + «2 + "3 )■ 2 + «2 + “3 )) " 2 ^2^p(" / + «2 + “3 ) 

(- (/ + «2 )+ 2 m2“^2 )■'■ ~2 “2 ^ ^3 ^1“^! ^ ^2 ^ 

d^ =-■ i-02/p(-2C2/ + / + 2M2 +M3)-yC202/p(-/ + M2 + Mj )~ y S2®2^p(^ + ”2 + ) 

(^66 =-3^2^2®2^^p(2“2 + 4«3 ) + -Lm2/p(- / + 2 C 2 / + 2U2 + Mj )+ y «3/p(- 2/ + 4CJ/ + 2 j + 2»3 ) 
^67= T®2^P(- 4/C2 + 2/ + «2 + 2U3 ) 

d - j-j = (/ + I/2)('y 02'Sl‘y2^P 2^ 02^1*'2^p)’^ T ^2^p{~ ^2^1 *^2)"^ ^*^2^1 ) 
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Eq. A14 


Q: 


Q = 


fxl ■*■ fx2 

fyl fy2 

'1 (- fxi^i + 4.^1 + “1 + “ 2 )' 

(-A 2 'Si+ 42 ^^iX^ + " 2 ) 

2 ” {fx\^l '*' fyl^l fx2^2 fy2^2 j 

fx2^\ “*■ fy2^\ 

2 '("“ fx\^\ fyl^l ) 

t(- 42 - 52 + 42 ^ 2 X^ + “ 2 +“ 3 ) 
^ {/x2^2 ■*■ fy2^2 ) 


Eq. A15 
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Dynamics of a flexible, flying tether: 


The dynamics of a flexible tether are evaluated through a continuous cable system model 
undergoing axial vibration. This continuous system model considers longitudinal vibration 
before, during and after a capture episode, and during and after a release episode. The model 
system equations are derived from Fig. 1 as follows. 



The force balance for a small element of cable, dx is given as, 

- £4 ^ + SEA/^ + 3W,<&)+/(x,/)ac=p 


ar' 


dx Eq. 1 


with u = u(x,t) the tether displacement at location jc and time /, E, A, and />the cable modulus, 
area and mass per unit length. Assuming uniform tether properties and dividing through by dx 
yields, 

Eq.2 

This displacement function represents the flexible-body motion of the tether, while xr and xel are 
rigid body position vectors representing rigid body motion of the tether. f(x,t) represents the 
forces per unit length acting on the tether. 

Next, a solution to the flexible tether dispacement function, u(x,i) is constructed as a sum of 
mode shapes and time dependent generalized coordinates. 




Eq. 3 


r=l 


with Ur(x) representing the modes shapes of the tether. These mode shapes are derived as 
solutions to the spatial differential equation. 


as. 


C/(x) = C, sin(/&:)+ Cj cos,{px') 


Eq. 4 
Eq. 5 



Eq. 6 


with B = — — . Coefficients C; and C 2 are determined from the boundary conditions to the 
EA 

problem. For the free-flying MXER tether, consider the homogenous boundary conditions 
associated with free-free tether ends. These are described as, 

£4«"('.')/=0 

Solving Eq. 5 for coefficients Cj and C2 yields the following mode shape vector and natural 

U, (x) = cos(^) Eq. 7 

Eq. 8 


frequencies. 




EA 


for r = 0, 1, 2, . . . with r = 0 representing the rigid body mode of the tether. 

The modes in Eq. 7 are normalized resulting in the normal modes for the free-free tether as 

Ur W = -j^cos(^), r = 1, 2, . . . 

with the following normalizing functions employed, 

{xp, {x}k = S„rr,s=\,2,-.- 


that also satisfy the function 




Eq. 9 

Eq. 10 
Eq. 11 


With the mode shapes derived, the displacement solution form of Eq. 3 is substituted into Eq. 2 
to result in, 

CO 00 w 

ffjx,. +waxx,, +7«o)XG)xx,, + + /wax5^t/^(x)g,(/)5 + wcD X o X ^ 17, (x)^, (r)2 

>-1 £q 12 

+ 2/mo) X (r)l + X 0)^ 

r-1 r=l ^ 

or 

(^)?r (^)^ + 2 W CO X t/, (x)^, (r)? + m(o X +CO X 1» x) J 17, (x)g, {t)u 

r=l r=l r=l £q 13 

+ YgXt)EA^-^^u = f(x,t)-wj(x, +OXX., +(0X<DXX,,) 

with a, Q the rigid body rotation rates and accelerations, x, a reference vector locating the tether 

(fig. 1) and Xei a vector running along the tether from the reference point to the location of the 

element of interest. Using the fact that the modes Ur(x) have been normalized based on Eq. 10 
and satisfy Eq. 11, Eq. 13 can be multiplied through by Ur(x) and integrated over the tether 
length, L as. 


f + 2 /w(d X ^ f + ^(o x +e> x eo (xX/,(x)?,(r>irM 

0 rt=l ^ ^0 

^ JV, (x)EA ^ {()^ = jpr {xX^M - + a X I,, + (0 X fi) X X,, ))& 


r=0 

+ 


Eq. 14 
Eq. 15 


To yield the following set of differential equations, 

{t)u + 2<oxq^ {t)u + +(a x +o) x © x)q^ {t)u + 0)^q, (t)0 = Q, {t)Q, 

This vector set of equations can be decomposed into axial and transverse components as, 

?,(r)+(-e=(0+®,k(') = a.-^(') Eq. 16 

the axial components and 

e(fM>)*0(‘M‘)=Q,.^(') Eq. 17 

with coh given in eq. 8, 0{t), &{t) the scalar tether rotation rate and acceleration respectively, and 
Qr oxiaM\Qr iroM (0 transverse external forces and rigid body inertial effects 

respectively with (f) given as, 

Qrif) = 0- + axx.,+©x®x x,,)>& Eq. 19 

forr= 0,1,2,... 

The above equations, for example Eqs. 15 and 19 demonstrate the inherent coupling between the 
rigid body motion and the elastic response. The response of the rigid body tether is considered 
next. 


Response of Rigid Body Tether: 

The rigid body tether motion under the effects of a spherical gravity model is considered next. A 
uniform, rigid tether is assumed in a planar orbit, as shown in fig. 2. 


Figure 2; Uniform rigid tether in planar orbit. 

The energy of the system with respect to the generalized coordinates is given as, 

V = ZEHl 

n 


Eq. 20 


with m,, the tether mass. It the tether mass moment of inertia, rt, Vt and ^ the position, velocity 
and angular rate of the tether center of mass. 


Eq. 21 


r, =x^e'* 


V, = x^e'* + x^i^'* +^/^' 

Lagrange’s equations for the rigid body tether result as. 


m,x. -m,x,6^ + 


fjm 






m.x, 0 + 2m,x.xd + 


2 /2 / 


^cos{0-(f^ 


= 0 


xjs\n{0 




= 0 


(«,'X+a)+ 


//wi,(x^/sin(i0 


(>^r' + cos(^ - 


= 0 . 


Eq. 22 
Eq. 23 
Eq. 24 


Equations 22-24 present the tether rigid body motion, while eqs. 3 and 16 give the elastic body 
displacements of the tether. A method of evaluating the total elastic behavior of the tether 
system is now considered. 

Consider a tether system as depicted in figs. 1 and 2 starting with initial orbit and rotation 
conditions, Xr=6000000e'", &o=0, thetadot0=.13, thetaddot0=0. First, the rigid body motion of the 
tether is evaluated, allowing the tether body angular rate to be observed. Equations 22-24 are 
evaluated over a period of time with solutions for Xr, 0 and ^ evaluated. Plots of q, 0{t\ 9(f) are 
shown in fig. 3 

tether angular rate 



Fig. 3, plot of example rigid body tether angular rate under steady-state conditions 

From fig. 3, the angular rate can be approximately described by the following sinusoidal 
fimction, 

9^ - c, cos(fi>r)+ Cj 


Eq. 25 




or as a taylor series expansion 


(-1)" 


9^ = VI 

*^(2« + l) 


+ c. 


Eq. 26 


with ci=. 0075, C2=.00275 and w = pi/75. 

The general elastic response for u(x, i) for a set of initial conditions can now be presented. First, 
consider the solution for the generalized coordinates of Eq. 16. The homogenous solution to Eq. 
16 is: 




^2 “ ^1 


1 - 


+(c, +a>jy 


+ ... 


X- 


C, +0, 


X' +- 


40 


(c, +(oJ 




Eq. 27 


with the general solution of q given in the form 


Eq. 28 


n=0 


and ao, aj coefficients determined from the initial value problem. Using these solutions for the 
homogenous problem and the method of variation of parameters, the particular solution for q is 
given as 

^ (^k'2 W - 9 1 {^)qi (^) ^ <ix (^)? 2 W- ^ . (^)?2 (^) 

Finally, the elastic displacement of the tether is given as. 










.Eq. 29 


The solution for u(x,t) depends on the generalized forcing function, Qr. This depends on both the 
external forces on the tether and the rigid body tether motion. Consider the distributed force to 
be a constant gravitational load, fo and the rigid body motion given in the following section. 
Then, the generalized forces become, 

QM)^ ••• 

or 

E,.20 

substituting the value of Qr back into eq. 17 yields. 




r=l 


9i(^V 2 W-^'ltk2W 


dr 


This equation can be expanded for the first relevant modes and evaluated 


